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Abstract 

We investigate the possibility for a fiat Bianchi I brene Universe to rec- 
ollaps due to the presence of a negative " dark radiation" and an anisotropie 
stress in the forni of homogeneous magnetic field, localized on the brane. 
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The cosmology of a brane Universe, intensively studied during last several 
years, has a variety of possible regimes. Some of these are impossible in the 
standard cosmological scenario. The quadratic stress-energy term in the right- 
hand side of the effective equation of motion ^ [2] can result in a nonstandard 
effective equation of state of a brane matter p — (7 — with 2 < 7 < 4, 
which is forbidden in the standard cosmology. Another important feature of 
brane models is the presence of projection of 5-dimensional Weil tensor onto 
the brane in the equations of motion. This projection can be decomposed into 
scalar, vector and tensor parts [3]. The scalar part, the only compatible with 
FRW metric on a brane, also called " dark radiation" . It can be either positive or 
negative, violating in the latter case the positive energy condition. Both these 
features (the quadratic corrections and the presence of the " dark radiation" ) can 
significantly modify the cosmological dynamics even in the simplest case of an 
FRW brane. One of the most important modifications alters future asymptotic, 



because a large enough negative " dark radiation" can cause a recollaps of a fìat 
brane Universe 0]. 

The dynamics of a Bianchi I brane also have been studied in detail. In 
general, as we lift the assumption of spatial isotropy, possible sources of an 
anisotropie stress should be taken into account. However, the tensor part of the 
Weil tensor projection onto a brane, being a nonlocal anisotropie stress, has no 
evolution equation [3]. When some form of nonlocal anysotropic stress is put 
into equation of motion " by hand" , the resulting dynamics significantly depends 
on the ansatz chosen|H|[7]. On the other hand, in the absence of the Weil tensor, 
a Bianchi I brane with a perfect fluid allows a complete description 

The next step in understanding details of brane cosmology have been done 
in [S], when a Bianchi I brane with a homogeneous magnetic field is studied. 
Being a source of anisotropie stress, the magnetic field is localized on the brane. 
The nonlocal anisotropie stress was chosen to be zero (so that the equations of 
motion can be written in a closed form), the "dark energy" 14, however, can 
not vanish identically, because the magnetic field enters as a source term into 
the evolution equation for U (see below). The resulting dynamical system has 
many stable points and a very complex behavior. Due to this complexity some 
important issues of magnetic brane evolution remain unclear even after a careful 
stable point analysis of 5 and require numerical studies. 

Denoting the combination of the gravitational Constant G and the brane 
tension A as a new parameter \ — (87rGA) _1 , one can write the set of equations 
which governs the evolution of a brane magnetic Universe in the form (see a 
detailed derivation and a description of the reference frame used in 0) 

e 2 = 3 P + + \h* + x fy + \ P h - h 4 ) + i8 X u (i) 
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h = --eh-2a+h (2) 

p = - 7 e P (3) 

= -ìe 2 -a 2 + ì[(3 7 -2)p+/, 2 ]-|[4(3 7 -l)p 2 + (6 7 +4)p/ l 2 +3^]-6 X W (4) 

Ù = ~&U-^[Bh 2 + \a + h 2 -^a +P \ (5) 

à+ = -0<7+ - a\ 2 - o\ 3 + \\?{\ - |[(3 7 - 2)p + (6) 

^- = -©^--^(^2-^3) (7) 

023 = — 0023 — 2(7i 2 CTi3 (8) 

Ó"12 = —0CT12 + (3(T+ + \/3cr_)(Ti2 + O13O23 (9) 

(713 = -0CT13 + (3(7+ - V3(7_)(7l3 + (712(723 (10) 



Here = 3iJ is the volume expansion rate, h is the nonvanishing component 
of the magnetic field (we use the reference frame in which two other components 
of the magnetic vector vanish), is the shear tensor. The variables + and 
a- are connected with the shear components as follows 

(CU, (7 2 2, O33) = (-2(7+, <7 + + V3(7_, (7+ - \/3(7_). 
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Moreover, it is possible to define new shear variables a a, <jb, &C using the 
transformation 



<Ji2 + icr 13 = a A e lv (11) 

V3a_ + ia 23 = (p B + ia c )e 2lv (12) 

Due to symmctries of our dynamical system, the equation of motion for the 
phase ip has the form 

cp = a c 

and the phase decouples from the system. The shear part of the dynamical 
system takes the form [5] 



&+ = -6a+ -a A + h 2 (l - |[(3 7 - 2)p + h 2 }) (13) 

ér A = (-6 + 3ct+ + a B )(T A (14) 

& B = -Octb -<? 2 a + 2e>c (15) 

&c = (-6 - 2a B )a c (16) 



Before describing our numerical results for the general system (l)-(5), (13)- 
(16) we summarize the conditions for recollaps in simpler models. 

If the matter content of the brane Universe is a perfect fluid only (so that 
h = 0) , the future recollaps conditions can be found analytically. First of ali, 
we can see from eq.(5) that U in this case behaves as radiation, U = Uo/a A , 
where a is a mean scale factor. The sign of ÌÀ remains unchanged during the 
cosmological evolution. 
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For an isotropie Universe the recollaps conditions have a very simple forni. 
Indced, for 7 > 4/3 the U - terni in eq.(l) falls less rapidly than the matter 
terms making the future recollaps inevitable for any initial negative U. 

For 7 = 4/3 the first and the last terms in the r.h.s. of (1) are equally 
important in a low-energy regime, and the condition for collapsc is p < —6xU. 
For 2/3 < 7 < 4/3 ali three nonzero terms in r.h.s. of (1) are important. 
Consider, for example, the case of 7 = 1. For the dust-like matter p = po/a 3 , 
and assuming that po, Uo are initial values of matter and "dark radiation" 
densities, taken when a = 1, we can substitut the scale factor dependences of 
p, p 2 and U into (1). As a result, we obtain that the condition H — lcads to 
a cubie equation for the scale factor 

3 , a M° 2 , XPO n 

a ó + 6x—a + -z- = 0. 
Po 2 

This equation has positive roots only for p^ < — 64% 2 Wq. A root which is 
bigger than 1 corresponds to the future of the Universe under consideration 
and indicates the point of recollaps. A root with a < 1 describes the past and 
corresponds to a bounce. 

A similar treatment can be performed for any other 7 > 2/3. If 7 < 2/3, 
even the quadratic stress-energy correction term, proportional to p 2 falls less 
rapidly than li, and the recollaps in the future becomes impossiblc. 

If we add a geometrical anisotropy without an anisotropie stress, the situ- 
ation does not change significantly. Indeed, the shear scalar S = 3cr^ + a 2 A + 
a B + a c behaves as E = Eo/a 6 , and we only have an additional term in the 
equation H(a) = 0. For example, the system with cosmological Constant A 
(which corresponds to 7 = 0) and shear gives 

So/a 6 + A + (x/2)A 2 + 6 X U /a 4 = 
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This gives us the condition for H = to bc possible as A 2 (l + (x/2)A) 2 < 
-32S x 3 W 3 . 

In a similar way we can treat other values of 7. As the shear E falls more 
rapidly than the dark radiation U, it becomes negligible in the low-energy limit, 
and for 7 > 4/3 with U < a recollaps remains inevitable. The condition for 
7 = 4/3 also remaines unchanged. 

If wc introduce an anisotropie stress, the shear and "dark radiation" can 
no longer be expressed in a simple forni as a function of the scale factor due 
to source terms in eqs.(5)-(6). In particular, now U does not proportional to 
a~ 4 and attribution of the name "dark radiation" to U becomes nothing but 
tradition. The problem of possible future recollaps of a brane Universe with a 
homogeneous magnetic field should be therefore investigated numerically. Our 
main goal in this paper is to describe qualitatively the influence of the magnetic 
field on the possibility of recollaps, keeping in mind that the non-compaetness of 
our phase space makes it impossiblc to introduce of a simple reasonable measure 
on it. 

Two main features distinguish the magnetic case from the case without 
anisotropie stress. 

• The /i 4 -term in (1) has minus sign, so the recollaps can occur even with 
nonnegative dark energy U if magnetic field h is large enough. 

• The sigh of dark energy U can change during the cosmological evolution. 

The sign change of U from positive to negative can lead to a recollaps for 
positive initial U even when h is not large (of course, in this case U < at 
the maximal expansion point). The opposite change from negative to positive 
values could in principal enlarge zones in initial condition space which do not 
lead to collapse. Both these possibilities are allowed by the equations of motion. 

However, in our numerical studies we have mostly noticed the transition from 
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Figure 1: Some trajectories which allow a transition from negative to positive 
values of " dark energy" 

plus to minus for the sign of hi. Some transitions minus — » plus were constructed 
using specially adjusted initial conditions, ali of them appeared to be part of a 
zigzag-like behavior and are not important for a general picture (see Fig.l). On 
the contrary, plus — » minus transitions are in some sense typical. As a result, 
in ali our simulations zones of initial conditions, escaping a recollaps shrink in 
comparison with similar models without a magnetic field. 

Some numerical results are shown in Fig. 2 for 7 = 4/3 using (if, /z) slices 
of initial condition space. In ali plots we put \ = 1- White zones represent 
the initial conditions leading to eternai expansion, gray zones lead to recollaps, 
while initial conditions from black zones are forbidden by the constraint (1). We 
can see that the zones, leading to recollaps enlarge and begin to contain some 
part of U > initial conditions as compared to H = 0. However, there is stili a 
part of initial state space providing an eternai expansion. 

Similar plots for 7 = 1 are shown in Fig. 3. 

Different plots, using the (h, /z) slices of the initial condition space with the 
initial U fixed are represented in Fig. 4. In particular, we can see how a recollaps 
being impossible for initial U > without a magnetic field becomes possiblc for 
nonzero h. 

As a result, we can see that a homogenios magnetic field, localized on a brane 
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Figure 2: The future fate of a brane magnetic Universe, starting the evolution 
from a a = 1-3, erg = °c = 1-1- Initial /i = O(lcft) and h = 2 (right). White zone 
corresponds to eternai expansion, gray zone - to recollaps, black zone represent 
initial condition, forbidden by the constraint equation (1). In this plot 7 = 4/3 

a 

The sanie as in Fig.2, but for 7 = 1 

diminishes the set of initial conditions, which results in a big Universe, resem- 
bling our oun. The set of initial conditions, leading to a recollaps, increases. 
However, this change is not drastic, and we can always find reasonable initial 
data for an eternally expanding Universe even if initial magnetic field is large. 
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